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It is known that in the regime of superlinear diffusion, characterized by zero integral friction 
(vanishing integral of the memory function), the generalized Langevin equation may have non- 
ergodic solutions that do not relax to equilibrium values. It is shown that the equation may have 
non-ergodic (non-stationary) solutions even if the integral of the memory function is finite and 
diffusion is normal. 
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There is hardly anything more important to say about a statistical mechanical system than whether it is ergodic or 
not. In general the question is notoriously difficult, yet for certain classes of stochastic systems the criteria of ergodicity 
breaking may be remarkably simple [lj45|. This is so, or so it would appear, for stochastic dynamics described by the 
generalized Langevin equation (GLE) 



dA(t) 
dt 



= - [ dr M(t - t) A(t) 
Jo 



F(t), (1) 



which governs a dynamical variable A of a classical system coupled to a thermal bath with many degrees of freedom 
in the absence of external forces [6[. The "random" force F(t) is zero centered (F(t)) = 0, not correlated with the 
initial value of A 

(A(O)F(t)) = 0, (2) 
and related with the dissipative memory function M(t) through the fluctuation-dissipation theorem 

(F(0)F(t)) = (A 2 )M(t). (3) 
We shall also assume the asymptotic vanishing of correlations of the random force 

lim (F(O)F(t)) = lim M(t) = 0, (4) 

t— >oo t— >oo 

which is typical for irreversible stochastic processes. It appears to be a common belief that, given conditions ([3]) and 
(j4|) , solutions of the GLE ([I]) describe ergodic relaxation to thermal equilibrium, unless the Laplace transform of the 
memory function M(s) = / °° dte~ st M(t) has a specific asymptotic behavior. Namely, it was shown in that the 
condition of ergodicity breaking for GLE systems has the form 

M(s)-s 5 , 8>1, as s^0. (5) 

This condition implies the vanishing integral of the memory function 

dtM(t) = M{0) = 0. (6) 



If the targeted variable A is the velocity of a Brownian particle, condition ([5]) corresponds to anomalous diffusion 
when the mean-square displacement (x 2 (t)) of the particle increases with time as t a with a > 1 (superdiffusion) 
The relation ([6]) is not very common, but not unrealistic. For instance, it was found to hold for a particle interacting 
with longitudinal phonons in liquids in the limit of zero temperature It should perhaps be noted that while the 
condition of ergodicity breaking |(5J| invariably implies superdiffusion, the converse is not true: for M(s) ~ s s with 
< S < 1 the condition of superdiffusion ([6]) is satisfied yet solutions of the GLE ([I]) are ergodic (0, see also Eq. (fT2|) 
below) . 

The purpose of this paper is to show that the condition of ergodicity breaking in the form ([5]) is too restrictive. It 
will be demonstrated that the GLE may have non-ergodic solutions even if the memory function does not follow the 
asymptotic form (O, M(0) = J Q M(t) dt is finite, and diffusion is normal. 
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We begin by briefly recapitulating the derivation of condition ([5]) which may differ depending on a type of averaging 
(...) in relations ^ and ©. When the GLE is derived with the Mori's projection operator technique @, the system 
is usually assumed to be in thermal equilibrium with the bath, and the averaging in equations (|2J| and (j3]) is over the 
ensemble of initial conditions for the composition of the system and the bath in mutual thermal equilibrium. In this 
case it is natural to use the GLE to evaluate the equilibrium correlation function (A(0)A(t)). Its normalized form 

_ (A(0)A(t)) 

c m- (A2(0)) (7) 



t 

dTM(t-T)C{r) (8) 



ii 



satisfies the equation 

dCit) _ 
dt ~ 

with the initial condition C(0) = 1, and has a Laplace transform 

C{s) = — (9) 
s + M(s) 

The connection to ergodic properties is given by Khinchin's theorem [tJ (see also Q), which states that the stationary 
process A(t) is ergodic if the correlation function factorises and, for a zero-centered process, vanishes in the long time 
limit 

fc^-W- 1 (10) 

Although in Mori's GLE the random force F(i), and therefore A(i), are not necessarily stationary and zero centered, 
we shall assume that these properties do hold. Then equations ([9]), (p~0|) and the limit value theorem 

lim C(t) = lim s C(s) (11) 

t— >oo s— >0 

give the condition of ergodicity breaking in the following form 

lim * ^ 0, (12) 

s^o s + M(s) 

which leads to condition ([5]). 

A slightly different approach is to apply for a particular and very popular class of models where the random force 
F(t) does not depend on A. This is the case, for instance, when A is the momentum of a Brownian particle which 
is bilinearly coupled to the bath comprised of harmonic oscillators [6]. For this problem, often referred to as the 
Caldeira-Leggett model, relation $2$ is satisfied trivially, the random force is stationary (for the infinite bath), and 
the fluctuation dissipation theorem takes the form 

(F(0)F(t)) = {A 2 )M(t), (13) 

where the the average (...)o is taken over bath variables only. The latter allows one to use the GLE to evaluate 
not only the equilibrium correlation function, but also the second moment (A 2 (t))o which characterizes the process 
of thcrmalization of the system which at the moment t = is put in contact with the equilibrium thermal bath. 
Compared to Mori's approach, this is a more general problem since the initial equilibrium of the system and the bath 
is not assumed. One can show that the system does not thermalize, 

hm(A 2 (t)) ^(A 2 ), (14) 

t— >oo 

under the same condition as that for ergodicity breaking discussed above. Indeed, using Laplace transformation the 
solution of the GLE ([TJ can be written in the form 

A(t)=A(fl)C(t)+ f dTC(t-T)F(r). (15) 

where the response function C(t) has the transform given by Eq.© and therefore coincides with the correlation 
function for Mori's GLE and satisfies Eq.©. By squaring and averaging solution (fl~5|) over bath variables, and also 
using stationarity of F(t) and the fluctuation-dissipation relation (fT3| 

(F(t 1 )F(t 2 )) = (A 2 )M(\t 1 -t 2 \), (16) 
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one gets 

(A 2 (t)) = A 2 (0) C 2 {t) + 2 (A 2 ) fdnCin) f Tl dr 2 C{t 2 )M{ Ti -t 2 ). (17) 

Jo Jo 

Using ((5J, this equation can be written as 



(A 2 (t)) Q = A 2 (0) C 2 (t) - 2 (A 2 ) f d n C( Tl ) C( n ), 

Jo 

and eventually one obtains 

(A 2 (t)) = A 2 (0)C 2 (t) + (A 2 )[l-C 2 (t)]. (18) 

The system does not thermalize if the response function does not vanish in the long time limit, lim t _ i . 0O C(t) 7^ 0, 
which again gives the conditions (fT2j) and ([5]). 

The above reasoning was based on the limit value theorem (|11[) which is only valid if the system reaches a stationary 
state and the long time limit for C(t) does exist. One might suggest that this is always the case provided the random 
force is irreversible in the sense that the correlation function (F(0)F(i)) and the memory kernel M(t) vanish as t — > oo. 
Let us show that this assumption is incorrect: It is possible to construct memory functions M(t) which vanish at long 
times, but the corresponding functions C(t), related to M(t) by Eq. (|SJ or ©, do not have a long time limit. The 
condition of superdiffusion J M(t) dt = is not required. 

As an example, let us consider a class of memory functions with the Laplace transform 

M(*) = -*+^p (19) 

where u is real, and f(s) is an analytic function at s = ±iu. The corresponding transform for the correlation or 
response function C(t), given by ((9]), is 

= 4tK- ( 2 °) 

It has simple poles at s — ±iui on the imaginary axis and therefore the original C(t) contains terms oscillating with 
frequency oj and does not reach a stationary value as t — > oo. It is not immediately obvious, however, whether it is 
possible to construct a function /(s) which ensures that the memory kernel behaves in a physically reasonable way. 
There are several conditions to satisfy. First, as a Laplace transform must vanish in the limit s — > oo, we must require, 
in view of (fT9"|) , that 

f(s) ~ s, as s — > oo. (21) 
The second condition is the asymptotic vanishing of correlations 

lim M(i) = lim s M(s) = 0, (22) 

t— >oo s— >0 

which leads to the asymptotic constraint 

f(s) - s r , r < 1 as s -> 0. (23) 

The third condition, which is more difficult to handle than the other two, is that the memory function must not 
exceed its initial value, 

\M(t)\<M(0), forOO. (24) 

This is because M(t) is essentially the correlation function of the stationary stochastic process F(t), which satisfies 
the inequality 

<[F(0) - F(t)} 2 ) = 2(F 2 (0)) - 2<F(0)F(t)> > 0. 

Condition ([24]) cannot be formulated as that for the Laplace transform M(s), which makes the choice of f(s) in 
Eq. (|19[) not quite straightforward. It turns out that the simplest function /(s) that can be made consistent with all 
three conditions $ZT$, and is 

/(-) = ^±», (25) 
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with certain restrictions on the real constants a, b, and c. In this case 



C{S) = {S + CW+U,*) ' 



while the transform of the the memory function (|19[) can be written in the form 



M(s)=a+ 2 P * + \, , (27) 
s + as + b 

where 

a — c — a, (3 = oj' 2 — b — a(c — a), 7 = ui 2 c — b(c — a). (28) 



The inverse Laplace transformation £ 1 of (|2T[) gives the memory function as a sum of the Dirac delta-function and 
a non-singular part, 

M(t) = a5(t) + m(t). (29) 
Conditions (|2"21 and are satisfied if the singular part is positive (a > 0) and the non-singular function 



m(t)=C-^ f S + 1 A (30) 
1 s 2 + as + b< 



vanishes as t — > 00. This sets the constraints 

c>a>0, b^O, (31) 

while uj is still an arbitrary parameter. 

As an illustration consider the set of parameters a = b = 1, c = 2, and = 2. Then equations (|28|l give a = 1, 
/3 = 0, 7 = 3, and transforms (f2l)|) and (|2"T|) read 

*<»> = 1 + :?Tbr 6 M= (/ + + 2X, + + V < 32 > 

Respectively, in this case the memory function is 

M(t) = S(t) + 2v / 3e- t/2 sin hj- t^j , (33) 

and satisfies both conditions (|22|) and (p4|) . while the correlation or response function 

C{t) = - e- 2t + i cos (V2<) (34) 

does not reach a long time limit. Observe that the non-singular part of the memory function m(t) increases at t = 0. 
As one can check, for the given class of memory functions (IT^l) . and under restriction Q, the property m'(0) > is 
generic. Therefore the presence of a singular term in M(t) is essential: if the delta function is absent in (|29|) (a = 0), 
then the condition (|2"4")| cannot be met. 

Needless to say, the condition of superdiffusion J °° dt M(t) = M(0) = 0, is not implied in our construction. It 
follows from (|2"?| that M(0) = uj 2 c/b, and so, unless w = 0, M(0) is finite and diffusion is normal. Since b ^ due to 
([3Tj). the regime of subdiffusion lim s _ ! .o M(s) = 00 Q does not occur for memory functions of type (|2"T1) . 

Summarizing, it is shown that stochastic dynamics governed by the generalized Langevin equation may be non- 
dissipative in the regime of normal diffusion, and thus superdiffusion is not a necessary condition for ergodicity 
breaking, as often assumed in literature. In our showcase example the memory function consists of a delta-peak and 
a long, generally nonmonotonic tail. A possibility of non-ergodic dynamics generated by a noise with physically more 
realistic non-singular correlations remains to be examined. In this case the time-reversal symmetry requires that the 
exact memory and autocorrelation functions must be even in time [l2l ]. Neither this additional constraint, nor the 
condition of subdiffusion (divergence of M(s) as s — > 0) can be met by the simple class of memory functions considered 
in this paper. 
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